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Abstract 



We consider a simple cosmological model that includes a long ekpyrotic contraction 
stage and smooth bounce after it. Ekpyrotic behavior is due to a scalar field with a 
negative exponential potential, whereas the Galileon field produces bounce. We give 
an analytical picture of how the bounce occurs within the weak gravity regime, and 
then perform numerical analysis to extend our results to a non-perturbative regime. 

1 Introduction 

The Universe possesses several features that are difficult to explain within the hot Big Bang 
theory. An elegant solution to these problems is inflation [U [21 EH E] ■ However, alongside 
with the inflationary theory there exist a number of cosmological scenarios which can also 
help to deal with the hot Big Bang theory problems. One set of ideas has to do with 



ekpyrotic El El El E3 , bouncing [TT1 US Efl H and cycling [151 ESI El US] cosmologies: 



long enough contraction epoch preceding the hot Big Bang stage could be an alternative to 



One way to consistently obtain the contracting pre-Big Bang stage is to introduce a scalar 
field with a negative exponential potential: V$ oc — exp(— af>). This kind of potential leads to 
the stiff equation of state, p/p > 1. Stiff matter avoids the "chaotic anisotropy" effect [19]; 
otherwise one would encounter anisotropic collapse of the Universe which is hard to convert 
into smooth expansion. On the contrary, contracting stage with stiff matter is consistent 
with isotropic, spatially flat and homogenious Universe [IS] . 

For the bounce to occur one has to ensure that the initially negative Hubble parameter 
crosses zero. Its derivative is given by the Raychaudhuri equation: 



inflation. 



H 



4irG(p + p) 
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so that the bounce requires 

p + p < 0. 

In other words, the bounce needs the dominating matter that violates the Null Energy 
Condition. The problem is that most types of matter do not behave in this way, while 
N.E.C. -violating theories are often plagued by istabilities. One known exception is Ghost 
Condensate, whose instabilities at the bounce stage may be mild [5J H3] and, indeed, a 
model with initially stiff equation of state and bounce due to Ghost Condensate has been 
recently constructed [H] . Another way is to introduce the Galileon field (2TJJ, EH [22] , or its 
close relative, a field with Kinetic Gravity Braiding [23l 124"] . These are scalar fields with 
special forms of the Lagrangian that lead to N.E.C. -violation but allow to avoid instabilities. 
Bouncing cosmological models with these fields have been proposed in Refs. [2S1 l2"5] . 

Our purpose is to construct a cosmological model which describes the Universe that starts 
from ekpyrotic contraction stage and proceeds to expansion via smooth bounce produced by 
the Galileon. In particular, we make sure that the cosmological constant is zero at late times. 
We find that ekpyrotic and bouncing behavior is quite generic, as it occurs in a wide range 
of the parameters and initial conditions. 

The paper is organized as follows. We start Section [2] with the description of a simplified 
model, then we study the system behavior at early times (Section 2.2) and consider the 
Galileon that evolves in the background of the scalar field. In Section 2J3 we introduce 
the perturbative approach and observe the bounce analytically, whereas in Section 2.4 we 



present numerical analysis. In Section [3] we modify the model to ensure that the cosmological 
constant vanishes at late times and show that the late-time inflationary behavior is avoided 
in a wide region of the parameter space. 



2 Bounce in a simplified model 
2.1 The Model 

We begin with a simple model with ekpyrosis and Galileon bounce; the model is modified 
and made more realistic in Section [3j We introduce two scalar fields, the Galileon 7r and the 
second scalar field 0, both minimally coupled to gravity (mostly negative signature): 

S = Seh + Sfj) + S n , 

where 

oeh = 7T~ / d X riy—g , 



S„ I <i l .r v'-fjilfTd^t-Vty) 
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Here (Ott 



\2n 



(5 ,A " / <9 /i 7n9 I ,7r) ri , n = V 2 . This particular version of the Galileon generically 



has a problem with super luminal propagation. However, this property is absent near the 
relevant backgrounds in a slightly modified model [22], in which the coefficient of the last 
term in @ is (1 + a)f 3 /2A 3 instead of f 3 /2A 3 , where < a < 3. To simplify formulas 
below, we set a = 0, but the analysis goes through for subluminal a. 
We study spatially homogeneous fields in FLRW background 

ds 2 = dt 2 - a 2 (t)dx 2 . 

Since we want to start from contracting Universe, we take the potential of (J) in the following 
form: 

V{<j>) = -V Q e~ c * , 

where Vo and c are positive parameters. Energy densities and pressures of the scalar field 
and the Galileon are 

P<t> = 7T + V{<f>) , 



Ptt 



fe 2 ^ 2 + 



f7T 4 + 4H7T 3 



f 2 

Pn = -f e 7T + 7^p{n 



Att 2 tt 



Note that the potential V((f) is unbounded from below. We improve on this point in Sec- 
tion [3l 



2.2 The Galileon in the background of scalar field 

If the Galileon field starts at 7r — > — oo as t — > -co, while the field evolves from <fi = oo, 
the Galileon energy density is sub dominant, and the cosmological contraction is governed 
by the scalar field. In this section we study ekpyrotic contracion of the Universe and the 
Galileon field in the contracting background. We are not able to obtain an exact solution 
for the Galileon; however, a combination of analytical and numerical analyses helps us to 
understand its behavior before the bounce. We will use these results in the next section as 
the initial conditions for the "slow bounce" scenario. 
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If one switches off the Galileon, it is straightforward to obtain the background solution 
for the scalar field and scale factor: 



<j)(t) = — In 
c 





2 1 


K 1 M p ¥y 


c 2 V t 2 _ 



a(t) oc (-t) c ' M p 
H(t) 2 



c 2 M 2 t 



To study the Galileon behavior, we write down its field equation: 

e 2w {n + tV 2 ) - -^-k 2 n + 3e 2w Hir - At UHtttt + 2ix 2 {3H 2 + H) + 2Htt 3 
H 2 H 2 L 

Here 



0. 



2A S 



In the case of Minkowski background, the stable solution is [21 

1 



H* (i* — t) 

It is parametrized by a parameter t*. In the case of evolving background with 

P 2 
H=-, P 



t 



c 2 M 2 



it appears natural to search for the solution in the form 



HA 



with some constant a. One finds that eq. ^ is indeed a solution to eq. Q with 



a 



6P 2 -8P + 2 
2-3P 



(3) 



(4) 



(5) 



However, this solution is a repulsor. To see this, we consider small perturbation n — > ty + 5ty. 
The linearized equation for 5n has the following form: 

t 2 5n(a 2 - 2 + 4P) + ^(4 - 2a 2 + 3a 2 P + 12P 2 - 2P) + 5it(Aa 2 - 6a 2 P) = . 

To simplify calculations we assume that c ^> M~ , work to the linear order in the small 
parameter P = 2(cM p )~ 2 and get the following solution: 



Sir 



C-i 



In the case of Minkowski background (M p — > oo, P — > 0), the growing part of perturbation 
is proportional to 1/t. It can therefore be viewed as the result of timeshift: if one makes a 
shift t — > t + e, the Galileon transforms into 7r — > n + 7re = 7r — e/t. In the case of curved 
background the growing part of perturbation cannot be interpreted in this way; this fact 
leads to the conclusion that ^ is a repulsor rather than attractor. 

The true solution in evolving background with P <C 1 behaves as follows. At early 
times the terms with the Hubble parameter are negligible in eq. Q, and the solution is the 
Minkowski attractor, 

i 

t -oo. (6) 



H*(h-t) 

Later on, the effect of the Hubble parameter accumulates, and the solution deviates from 
eq. ([6]). However, after a certain transition period, the Galileon field and its time derivatives 
become so large that the Hubble terms in eq. Q become negligible again. The solution 
asymptotes to 

-£/* (t* — t) ^ 
with some new parameter t*. The latter regime occurs when 

\U-t\ < ^. 

An example of numerical solution is shown in Fig.l. 



" H(t) 




Figure 1. The Galileon in the background of scalar field. Numerical solution (solid line) and 
asymptotics e n = —^~ t (dashed line) and e n = H } _ t , (dot-dashed line). 



The solution ([7]) is valid modulo small corrections, provided that the terms with the 
Hubble parameter are small in the Galileon field equation. As we will now see, this situation 
can occur even at the bounce epoch. We call this regime "slow bounce". 

2.3 Slow bounce 

The Galileon equation has rather complicated form and the full system of equations cannot 



be solved exactly. However, there exists a useful approximation. We have seen in Section [2^2 
that at times of interest the Hubble terms in the Galileon equation are negligible, and the 
Galileon rolls just like in Minkowski space-time, see eq. ([7]). Likewise, for M p c 1, the 
scalar field (|3| also exhibits Minkowskian behavior. This suggests that we can make use of 
perturbative approach in the small parameter M~ 2 . 

To the leading order in M~ 2 , the Galileon and scalar field <fi decouple. We use their 
Minkowski space solutions, eq. ([7]) and 

2 



c c 2 V t 2 ' 

and find that energy densities of the Galileon and scalar field vanish, while their pressures 
do not: 

PL 0) = 0, P^O, (9) 

^ o) =-3 (£)%*-*r 4 , 



We then use eq. Q to calculate the first-order Hubble parameter: 

^=M* + m*±v)- (10) 

Here the integration constant is set equal to zero because of the initial condition H — > as 
t —7- — oo. 

Note that H^ 2 oc M~ 4 , and hence, in view of the Friedmann equation, pc + oc M~ 2 . 
This is consistent with eq. ([9]): energy densities become non- vanishing at the first order of 
the expansion in M~ 2 . 



As it stands, eq. (10) describes ekpyrotic contraction and bounce. We have to check, 



however, that bounce occurs within the range of validity of our perturbative treatment. 
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2.3.1 Conditions of validity 



Our perturbative approach is valid provided that the gravitational corrections to the field 
equations are small. We are interested in the regime when H crosses zero; that happens at 
time tb such that 

2 p 1 



+ 



0. 
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is 



cH b 3#2 (U - t b ) 3 

It is straightforward to see that the gravitational correction to the field equation for 
small provided that 

c 2 M 2 p > 1. 

The terms in the Galileon equation Q that do not include gravity are of the order [H 2 (£* — tb) 
at t = tb, whereas gravity corrections make two combinations: 

[Hlc%{t h -UfY\ [HlcHHtb-UfY 1 - 



(12) 



In view of eq. (12) they are small, provided that 

cM p \t b \ > \U-t b \. 



(13) 



Let us now study solutions to eq. (11) to see whether eq. (13) is satisfied. Note that our 
solution has two singular points, t = t* and t = 0, so all physics occurs at t < 0, t < t*. 



First, one can check that for t* < 0, eq. (11) has one solution in the t < region, and this 



solution satisfies the condition (13), so the bounce inevitably occurs. For £* > 0, there are 
two possibilities: 



1. £* ^> cf/H*. In this case eq. (11) does not have negative solutions. This means that 



there is no bounce within the range of validity of our perturbative treatment. 

2. t* <C cf/H*. This inequality implies that there exist two negative roots of eq. (11) 



One of them is always within the region of validity and corresponds to bounce, whereas the 
second one describes "rebounce", after which the Universe contracts again. The rebounce 
occurs at 



c 2 f 

This second root satisfies eq. (13) provided that 

P 



2 L * 



tl> 



(14) 



If the relation (14) is violated, there remains one root of eq. (11 ) within the region of validity. 



and simple bounce occurs. 

Three possible types of behavior we have found within our perturbative approach are 
shown in Fig. 2. 
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2.4 Numerical results 

We have studied numerically the full system of equations describing our model. In the 
first place, we have checked that numerical solutions coincide with the analytical results of 
Section 2.3 within the range of validity of the latter. In particular, once simple bounce (no 
rebounce) occurs when the perturbative treatment is applicable, the Universe continues to 
expand in the non-perturbative regime as well. This is illustrated in Fig. 3. 



numerics 




Figure 3. Analytical solution (solid) vs numerics (dotted) 



However, numerics does not show two other types of behavior: "no bounce" and "re- 
bounce" with subsequent collapse. The bounce is always there. To see what is going on, 
let us choose some value of t* > in eq. ^ treated as the initial condition, and vary the 



parameter /. For small / (case 1 of Section [2^3 ) there is no bounce analytically, but the 
bounce does occur outside the range of validity of the perturbative treatment (Fig. 4a). For 
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larger / (case 2 of Section 2.3 with eq. ( [14] ) satisfied), there is indeed a rebounce, but after 
the second contraction stage the Universe bounces again, now in the non-perturbative regime 
(Fig. 4b). For even larger /, when eq. (14) is violated, the rebounce does not occur and the 
Universe continues to expand after bounce (Figs. 4c, 4d). 



H(t) 





H(t) 



(a) (b) (c) (d) 

Figure 4. Numerical results. Behavior of the Hubble parameter depends on the Galileon 
strength that is determined by the parameter f 

At qualitative level, the fact that the Hubble parameter does not decrease indefinitely but 
starts to grow instead, can be understood in the following way. Suppose that H{t) decreases 
and has large negative value, so that it dominates in the Galileon equation One observes 
that to the leading order in H, the Galileon equation reduces to 



7T = — t\H. 
2 

In this approximation, the energy density and pressure of the Galileon field are 



(15) 



p 7: . I - , / ' /2 "''// 



We plug these expressions into eq. Q together with (p + p)^ = 2 and obtain 

H = -AixG (q^tt 3 H + 2 

Since the Hubble parameter has large negative value, the system eventually enters the regime 
in which 



< |6^-7i- 3 iJ|. 
Hi 



(16) 



In this case the derivative of the Hubble parameter becomes positive and H begins to grow. 
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To substantiate this qualitative argument, let us check that the regime (16) is self- 
consistent. To this end we take the derivative of that inequality: 



f 2 

200 < \ 4-(37r 2 nH + tt 3 H) 
rii 



(17) 



Then we use the field equation for together with eq. (15), and find that eq. (17) transforms 



back into eq. (16). We see that eq. (16) is consistent with field equations: its left hand side 



grows slower than the right hand side. Hence, once our system has entered the regime (16), 
it does not leave it. 

To end up this section, we note that in a wide range of parameters and initial data, all 
regimes we have found take place at sub-Planckian energy densities and pressures. However, 
pushing the parameters and/or initial data, one does find situations when the bounce does 
not occur before the system reaches the Planckian parameters. This happens, in particular, 
for very large positive values of the parameter £„ characterizing the initial data for the 
Galileon: at the level of classical field equations, the bounce does occur in this case but it 
cannot be trusted, since pressure exceeds before the bounce. 



3 Extended model 

In the — 7r system, the Galileon dominates at late times. To be more realistic, one has 
to assume that the Galileon eventually decays, i.e., its energy density gets transformed into 
heat, see Ref. [27J • We also have to improve on the dynamics of 0, since, as we mentioned 
above, the potential V(<j>) is unbounded from below. We modify this potential in such a 
way that it obtains a minimum at finite 0. Our analysis in Sec. [2] will be valid, provided 
that the potential remains nearly exponential before and at the bounce epoch, so that our 
modification affects late time dynamics only. 

Let us denote the depth of the modified potential by V min ((f>) = —Wo- To make sure 
that the cosmological constant vanishes at late times, we also introduce a potential for the 
Galileon, W(ir), that asymptotes to Wo- The picture of Sec. [2] holds if it is equal to zero at 
the stage of bounce and earlier, but grows to the value Wo at late times. In fact, we have 
seen by numerical analysis that the whole picture often does not change even if the Galileon 
reaches W before the bounce. 

In the modified model, the field does not indefinetly roll down. There are several 
possible types of its late-time evolution. The first one is that reaches the minimum of 
V(4>) and settles down there after oscillating for some time, while the Galileon climbs its 
potential and decays later. The second possibility is that due to the rapid growth of the 
Hubble parameter (which is driven by the rapidly growing Galileon field), the scalar gets 
frozen before it reaches the minimum of V(4>). In that case its evolution stops at some value 
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4> s and continues only after the decay of the Galileon field. This situation can in principle 
lead to the standard inflationary stage, but we will show that this does not happen in a wide 
range of parameters. 

Let us discuss the dynamics using a specific example of modified potentials. 

3.1 Example 

Let us choose the scalar potential in the following way (left graph in Fig. 5): 

V(<f>) = -Voe-^ [1 - ae~ H ] , b > c. 
We also introduce the potential of the Galileon field (right graph in Fig. 5): 



W(ir) = 




[1 + tanh(g(vr - m))} . 



Here q ~ 1, and 7Ti is the parameter that determines the time when the Galileon starts to 
climb the potential. We choose W to be equal to the depth of the modified potential V((f)), 
so the cosmological constant vanishes at late times, 




-> 7T 



Figure 5. Modified potentials and the original potential V(<f)) (dotted line). 

Numerical analysis shows that depending on the parameters, the scalar field <fr can either 
reach the minimum of its potential V(<f)) and oscillate for some time or stop earlier due to 
the Hubble friction. The first situation does not require further analysis: near the minimum 
of V(4>) we are dealing with massive scalar field whose value is well below M p ; hence, the 
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model never enters the inflationary regime. We will shortly give an example of the parameter 



choice yielding this situation, see eq. (22). 



Less trivial is the case in which the field does not reach the minimum of its potential 
and becomes frozen at some point <p s = <f>(t a ). Let us see that in a wide range of parameters 
the inflationary slow roll conditions for V(<f>) are violated, so the inflation stage does not 
occur after the decay of the Galileon. To this end, we consider initial conditions such that 
t* ^> / /M p H* for the Galileon, and obtain an estimate for t s . The scalar field equation is 

<j> + cV e~^ + 3H<j) = 0. 

The field <ft begins to slow down when 

101 ~ |3#0|- 



Plugging the Minkowski solution (|8j) for the field (f> into this expression, one gets for freeze-out 
time 

H . (18) 

ts 

This occurs when the Galileon dominates in the Hubble parameter (we check this below), so 
that 

f 2 1 

H ~ — (19) 



Combining eqs. ( 18 ) and ( 19 ) one gets 

P 



M*H*\U-t a \* \t 



(20) 



Under our assumption £* ^> f/M p H*, we find that eq. (20) is satisfied when |t* — t a \ <C £*, 
so that 

Let us check that at this moment the slow bounce regime is still valid for the Galileon 
field. As we pointed out in Sec. |2.3.1[ the terms in the Galileon equation that do not 
include gravity are of the order [H^(t* — t) 4 ] . Corrections due to the Hubble parameter 
are H ■ [Hl(U — t) 3 ] 1 and H ■ [H*(U ~ t) 2 }' 1 . They are small provided that 

1 MlEl 

< p 



This inequality is indeed valid for t* ^> f/M p H*. 
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To check the consistency of our calculation, we observe that the scalar contribution to 
the Hubble parameter is small at t ~ t s . Indeed, we have 



> 



H*(U-t s ) 3 t s cH 



for M p c ^> 1, so the first term in the right hand side of eq. (10) is small compared to the 
second term. 

Now we turn to the main issue of whether the Universe enters the inflationary regime 
at late times after the Galileon decay. This can happen if the field freezes out at the 
exponential part of its potential, aexp(— bcj)) <C 1. This gives 

a(c 2 V t 2 s )-° < 1. (21) 

In the opposite case, 

t,< a 3F-L=, (22) 

the field stops very near the minimum of its potential, after oscillating about this minimum. 
As we have already pointed out, this behavior does not result in the inflationary stage. 
We proceed with the case (21). The inflationary slow roll parameter at = <p s is 



T) = M 



, V" (M p c) 2 (M p c) 



P V + W a-^e^ a- k b c 2 V t 2 s 
Therefore, the inflationary slow-roll condition is violated provided that 

a- £ bc 2 V t 2 s < (M p c) 2 . 

We recall that t s f» and conclude that the late time evolution of our system will not lead 
to the slow roll inflation at least for the following choice of initial data: 

/ « M p 



M P H* ^/Vo 

Clearly, if the parameters of the action for 7r and are far from Planckian, this range is 
parametrically wide. 



4 Conclusion 

We have shown that there exists a simple way to construct a bouncing cosmological model 
which incorporates a contraction epoch with p 3> p and bounce due to the Galileon. Notably, 
all this can happen in the regime where gravity can be treated as perturbation. This is not 
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a necessity, though: we have seen by numerical analysis that the Galileon dynamics is so 
strong that it always yields the bounce, irrespectively of the value of the Hubble parameter 
late at the contraction stage. We have made the model more realistic by ensuring that the 
cosmological constant is negligible after the putative decay of the Galileon and checked that 
the late-time inflationary epoch is avoided in a wide range of initial data. All this makes 
models of this type consistent, albeit somewhat baroque, alternative to inflation. 
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